Self-organization in ^He near the superfluid transition in heat flow and gravity 
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We investigate the nonlinear dynamics of ^He slightly below the superfluid transition by integrat- 
ing model F equations in three dimensions. When a superfluid is heated from above under gravity, 
a vortex tangle and sheetlike phase slips both appear near the bottom plate. Then a self-organized 
superfluid containing high-density vortices and phase slips grows upward, where high-amplitude sec- 
ond sounds are emitted from the self-organized to the ordinary superfluid region. A phase slip sheet 
often changes into a vortex aggregate and vice versa. The thermal resistance due to these defects 
produces a constant temperature gradient equal to the gradient of the pressure-dependent transition 
temperature T\(p). In this self-organized region, the temperature deviation T — T\{jp) consists of a 
negative constant determined by the heat flux Q and time-dependent fluctuations. Its time-average 
is calculated to be 155 nK for Q = 11.2 erg/cm'^s in good agreement with the experiment (W.A. 
Moeur et al, Phys. Rev. Lett. 78, 2421 (1997)). 

PACS numbers: 67.25. dj, 67.25.dk, 05.70. Jk,05. 70. Ln 



I 

o 



(N 

> 

o 

^' 

o 
o 



X 



I. INTRODUCTION 

The critical phenomena in ^He near the superfluid 
transition have been studied with exceedingly high preci- 
sion of the temperature T (~ 10~^ deg)i"— . Such experi- 
ments have thus been used to confirm or predict universal 
relations near the critical point. Moreover, in this sys- 
tem, a number of unique nonequilibrium situations can 
be found in the nonlinear response regime. In particu- 
lar, the superfluid transition is very sensitive to applied 
heat flow Qr^. For example, if T = T{z) is above the 
transition temperature Tx at one end of the cell and be- 
low T\ at the other end, a Hel-Hell interface emerges 
separating superfluid and normal fluid regional"—. Here 
T is nearly constant in the superfluid and has a finite 
gradient in the normal fluid, so the temperature gradient 
is almost discontinuous across the interface. In experi- 
ments, it emerges when "^He in a normal fluid is cooled 
from a boundary below T\ or when ^He in a superfluid 
is warmed from a boundary above T\. This interface 
has a thickness proportional to Q~^^^ at small Q and is 
thus a unique singular object. Some authors were also in- 
teresred in the heat capacity in heat flow near the super- 
fluid transitior>iS. Another interesting case is to quench 
liquid "^He from its normal to its superfluid phase by a 
mechanical expansioni^"— . In such phase ordering with a 
multi-component order parameter, topological singulari- 
ties (vortices in ''He) should be proliferated in an early 
stage after quenching^ii^. 

The superfluid transition can also be greatly influ- 
enced by gravity due to the pressure dependence of 
T\{p). In ''He near the transition, the specific heat ex- 
hibits a sizable logarithmic singularity, but the isother- 
mal compressibility has almost no singularity'^''^, result- 
ing in nearly homogeneous density profiles for not large 
cell length on earth. Thus the height-dependence of T\ (p) 
is the main origin of inhomogeneity in equilibrium *He, 
preventing precise measurements of the critical phenom- 
ena in ^He^''^. If the z axis is taken in the upward direc- 



tion on earth, T\{p) decreases with the height z as 

Ta(p) = rAbot(l + G2), (1) 

where TAbot is the transition temperature at the cell bot- 
tom and G = \dTx/dp\pg/Tx = 0.6 x IQ-^cm-^ at satu- 
rated vapor pressure with g being the gravity constant. 
The origin of the z axis is at the cell bottom. We intro- 
duce the local reduced temperature as 



T = T/Txip) - 1 = (T/TAbot - 1) - Gz, 



(2) 



which depends on z and represents the distance from the 
lambda line in the p-T phase diagram Therefore, equilib- 
rium states are noticeably inhomogeneous in the temper- 
ature region 



Irl < Gh, 



(3) 



where h is the vertical cell length. Remarkably, in the 
temperature range T{h) < < r(0) at homogeneous T 
(in equilibrium), a gravity- induced two-phase coexistence 
is possible with an upper region in a superfluid and a 
lower region in a normal flui d"' ''''^ . The thickness and 
the characteristic reduced temperature Tg of the gravity- 
induced interface is determined by the balance r = G^, 
where ^ = ■Cokl"^^'^ is the correlation length with = 
I.4A. Hence they are given by 



^9=e+o(e+oG)-2/5, Tg = (e+oG)3/5, 



(4) 



which are about 10^^ and 10~^ cm on earth, respectively. 
It is worth noting that application of an electric fleld 
E or a, magnetic field H gives rise to a shift of T\ by 
an amount proportional to E'^ or H^. If such a field is 
inhomogeneous, it leads to a gradient of Tx, as well as 
the earth gravityi§.. 

Furthermore, both gravity g and heat flux Q can be 
crucial near the lambda line. Hereafter Q will be mea- 
sured in units of erg/cm^s. If heated from below, they 
both serve to decrease r with increasing the height and 



2 



the interface changes over from the gravity-induced one 
to the heat-induced one at Q ~ 1 with increasing Q. (We 
may determine this crossover Q by equating the thick- 
nesses of these interfaces.) On the other hand, if heated 
from above, they can compete to produce new nonequi- 
hbrium states. In fact, Moeur ei al^— reahzed self- 
organized superfluid in ^He heated from above, where 
the temperature gradient and the gradient of T\ (j>) were 
balanced as 



az az 



(5) 



yielding a homogeneous t = — Tsos(Q) with 

r,os(Q) = 1.0 X 10-«Q°-8" (6) 

in the range 5 < Q < 60. This self-organized region 
extended over a macroscopic region in a 7.4 mm tall, 2 cm 
diameter cylindrical cell. It is surprising that the distance 
T to the lambda line was held fixed at an extremely small 
negative constant over a macroscopic region. Machta et 
al^ pointed out that ''He in heat flow is an example of 
self-organized criticality. However, ^He in heat flow and 
gravity stays slightly away from the lambda line, so it is 
not in a "critical state", to be precise. 

Notice that the balance relation (5) itself follows from 
the vortex resistance in the Gorter-Mellink forrc^, 



dT/dz = TB{T)Q^, 



(7) 



where B = So|t|-'^ with Bo = 5 x lO^^g .^^^ ^ ^ 2.23 
near the transitional. The balance B{T)Q^ = G then 
yields r 



-Cs with 



I/O 



1.3 X lO-'^Q 



10^1.35 



(8) 



However, this t^^j, is smaller than the experimental Tgos 
in Eq.(6) by one or two orders of magnitude for (5^5. 
In previous simulations in one dimensional^, another 
kind of singular objects, phase slips, came into play, pro- 
ducing a thermal resistance in accord with the experi- 
ment. In the literature, phase slips were theoretically 
studied in one dimension22r— and have been observed in 
quasi-one-dimensional superconductors^ and in ^He in 
narrow apertures^. The self-organized superfluid states 
are analogous to the resistive states in wire superconduc- 
tors in electric field; however, in our case, phase slips 
appear as sheets even for macroscopic cell width. This is 
because the uniaxial effect of gravity is intensified near 
the lambda line. Similar self-organization should well be 
realized in "^He in the presence of a gradient of electric or 
magnetic field, for example, around the tip of a charged 
needle kept cooler than the surrounding liquid *He. 

Originally, the balance relation (5) was predicted for 
normal liquid *He^. Liquid ^He exhibits strong criti- 
cal fiuctuations of the complex order parameter slightly 
above the lambda line, which give rise to strong critical 
enhancement of the thermal conductivity Ai^~. For small 



positive T, experimental data of A may be fitted to the 
following power law^, 



A = Aor 



{xx = 0.45), 



(9) 



where Aq is a constant. In terms of t in Eq.(2), the heat 
conduction equation reads 



dr 

dz 



G 



(10) 



where the specific heat Cp is treated to be a constant for 
simplicity. In steady states, the temperature gradient is 
Q/Xso that 



dr/dz = {Q/TxXo)i 



-G 



(11) 



If r(/i) > 0, the balance in Eq.(4) gives t = Tson(Q) with 

Tson(Q) = (TaAqG/Q)!/^^ = 3.8 X 10^90-2.22 (^3) 

This holds in the region L — z in{Q), where in{Q) = 
Tson{Q)/G is the relaxation length. Weichman and 
Miller— argued that the balance relation (5) is attained 
only when the correlation length ^ = Cokl^^^'^ is shorter 
than £n{Q)- From Eq.(4) this local equilibrium condi- 
tion is rewritten as Tgon <^ ''g ^ IQ-^. Thus we obtain 
the upper bound Q ^ 1 needed for the self-organization 
above the lambda line. For larger Q, the system crosses 
across the lambda line and a HeTHell interface appears. 
In accord with these results, Moeur et al}^-: observed self- 
organized normal fluid states for Q < 1. Analyzing the 
nonlinear diffusive equation (10), Weichman and Miller— 
also predicted a new thermal wave propagating only up- 
ward in self-organized normal states. This wave was ob- 
served both above and below the lambda line^^'^^. 

The physical processes in the self-organized super- 
fluid remain quite unclear, since in the experiments near 
the transition measured was only the temperature and 
the simulations were one-dimensional. Realistic three- 
dimensional simulations are needed to investigate how 
vortices and phase slips come into play. With this pur- 
pose, we will present dynamic equations in Sec. II and 
numerical results in three dimensions in Sec. HI. 



II. MODEL EQUATIONS 



^4.13 



We use the renormalized model F equation o"'"" ac- 
counting for gravity^. The dynamic variables are the 
complex order parameter il){r,t) and the entropy s{r,t). 
The superfluid density ps and the superfluid momen- 
tum density Js = Ps^s are proportional to li/jp and 
Im(V'*VV'), respectively, where Im(- • •) denotes taking 
the imaginary part. The Fourier components of ip and s 
have wave numbers smaller than a upper cut-off A. The 
fluctuations with wave numbers larger than A have al- 
ready been coarse-grained^ii^. As A is decreased down 
to the inverse correlation length the renormalization 
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effect is suppressed to justify use of tfie mean-field tlieory 
at long wavelengths, where the coefficients of the model 
are renormalized ones dependent on r. 

In equilibrium superfluid, the average order param- 
eter behaves as ■0 oc |t|^/'^. Without gravity, we de- 
fine the correlation length as ^ = ^+o| kl"^^^ using 
C+o = 1.4A. We suppose a reference equilibrium super- 
fluid state, where r = — f, ip = 'ip, and £, = £, with 



2.5 X 10" 



C = 1.6 X 10~'^cm. 



(13) 



Hereafter we set A = and treat the fluctuations with 
sizes longer than ^. Those with sizes shorter than ^ are 
little affected by heat flow. The scaled position r/^, the 
scaled time t/i, the scaled correlation length ^/^, and the 
scaled reduced temperature e/e will be written as r, t, ^, 
and e in the same notation for simplicity, where 



i = 0.5 X lO^^s 



(14) 



is the thermal relaxation time on the scale of ^. We 
introduce the scaled reduced temperatures. 



A{r, t) = {T/Txbot - 1)/t, 
e{r,t) = T/f ^A-Gz, 



(15) 
(16) 



where G = G^/f — 0.04 on earth. The equations for the 
scaled variables 5'(r,i) — ^jj/ij] and A{r,t) read^ 



-1/2 



r) A 

at a L 
d 

—M = alxa{^*VH) + V • A*Vyl, 



(17) 
(18) 



where M is the scaled entropy deviation defined by 

M = A~a^C^/^\^^l2 (19) 

and a is a universal number of order unity. The second 
term (cx l^'P) represents a decrease of the entropy in the 
ordered phase. We set a = 1 in our simulation. In gravity 
the correlation length does not exceeds £g in Eq.(4), so 
we define the local correlation length as 



^ = 4tanh(l/4|£|2/3) 



(20) 



where 4 = tg/i, = l/G^/^ ^ 3.62. Then we obtain two 
cases of C = \e\-'^/^ < Ig and ^ ^ < |e|-2/3. The 
scaled kinetic coefficients are written a.s^^ 



0.325 



A* = 6a^' 



0.675 



(21) 



We set 6^ = 0.2 and 6a = 1 in our simulation^. 

For "^He near the lambda line, the complete hydrody- 
namic equations including the dynamic equation for tp 
were presented by the Russian group many years ago23.. 
In the model F— , we do not treat the hydrodynamic equa- 
tion for the mass density p and the total momentum den- 
sity J = pnVn + Js neglecting the fluctuations of p and J, 
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FIG. 1. (color online). Top: scaled temperature deviation 
A{r,t) a.t X — y — after application of heat flux from the 
top at t = 700, 1000, and 1300 -I- 20n (0 < n < 4). The bulk 
region is gradually heated by repeated traversals of second 
sounds. Bottom: |\l/p(cx ps) (left) and d9/dz{(x Vsz) (right) 
a.t t = 900, 940, and 1340. A phase slip is being created at 
z ~ 10 in the latter two times (see the plates of Video 1). 



where Vn is the normal fluid velocity and p^ = p— p^ = p. 
Then the normal fluid velocity t>„ may be written as 



-p-'js 



(22) 



Therefore is much smaller than Vg = pj^-is by the 
small factor Ps/p- This is justified for slow dynamic pro- 
cesses without first sounds near the lambda line. 

In critical dynamics, the dynamic equations have been 
treated as Langevin equations with random source terms 
added*. In this work, we neglect the random source terms 
in Eqs.(17) and (18), which produce the thermal fluctu- 
ations of the order parameter -0 and the entropy s. They 
may be neglected when the small-scale thermal fluctua- 
tions shorter than ^ have been coarse-grained in the dy- 
namic equations. The renormalization effect has already 
be accounted for in the renormalized coefficients. 



III. SIMULATION RESULTS 

In previous papers^^i^, Eqs.(17) and (18) were inte- 
grated in one dimension, where phase slips play a major 
role in self-organized superfluid states, as in wire super- 
conductors in electric field'^^""^^. In this work, we present 
numerical results in three dimensions on a 200 x 200 x 200 



4 




■08 1,8 



FIG. 2. Surfaces of = 0.8 after application of heat 

flux from the top in a superfluid state in the same run as 
in Fig.l. Their colors represent the temperature deviation A 
with the color map below the plates. Top: a,t t = 940 a layer 
in the range 10 < 2: < 20 is much disturbed, where a phase 
slip is going to appear. Temperature inhomogeneity is also 
apparent, for which see Fig. 3. Middle: ai t = 1340 a sheetlike 
phase slip is created near the bottom. Vortices above the 
sheet are expanding. Bottom: at t = 2320 the phase slip 
sheet is away from the bottom and is being broken into an 
aggregate of vortices. Video gives time evolution in the range 
2000 < t < 3400. 



cubic lattice with the mesh size equal to 1^, so the corre- 
sponding cell length is 3.2 mm. The time increment At 
of integration is 0.01 in units of < = 0.5 x lO^'^s. The pe- 
riodic boundary conditions are imposed in the xy plane. 
At 2: = and 200, we set = and assume homogene- 
ity of the temperature and the heat flux in the xy plane, 
supposing metallic plates. We prepared an equilibrium 
superfluid state with A = —2 for t < 0. We then applied 
a constant heat flux at z = h, but we kept the bottom 
temperature at ^ = —2 (t = — 2f) at z = for t > 0. 
The scaled heat flux Q = X* (dA/dz) at the top is equal 
to 1. This corresponds to a heat flux of Q = 11.2, which 
is the value of XT\t/^ at t = f. In the following we will 



explain our numerical results. 

In Fig.l, we display profiles of the scaled reduced tem- 
perature A YS z (with X — y — 0) at seven times 
t = 700, 1000, and 1260 -t- 40?! (0 < n < 4), where traver- 
sals of second sounds cause gradual heating of the bulk 
superfluidii. Since the bottom temperature is fixed, a 
transition layer supporting a temperature gradient grows 
near the bottom, where a heat flux of order 0.1 is from 
the liquid to the bottom plate. This heat output is 10% 
of the heat input from the top. In the lower plates, we 
show l^fp and d9/dz near the bottom at three times to 
illustrate formation of a phase slip, where 9 is the phase 
of vfi = l^ile*^. We confirmed that this phase slip forma- 
tion in a narrow region occurs even if the mesh size is 
reduced to 0.2^. Here dO/dz is equal to the superfluid 
velocity in the z axis in units of H/m^,, where m is the "^He 
mass. Note that the complex order parameter behaves as 
■0 oc e''^^ with Vg — hk/m in current-carrying superfluids 
and the bulk region is unstable for k larger than 1 / \/3^ 
from linear stability analysis^. In our case, ps decreases 
towards zero and Vs exceeds the critical value near the 
bottom, resulting in defect formation. 

In Video 1, we display the surfaces of — 0.8 

near the bottom, which illustrates growth of vortices and 
phase slips. At t — 940 disturbances of 5* are enhanced 
in a layer in the range 10 < 2; < 20. See the lower 
plates in Fig.l. We can also see considerable inhomo- 
geneity of the temperature. At t = 1340, vortices are 
expanding and a sheetlike phase slip is at its birth. In 
our case under gravity, Ivfp decreases near the bottom. 
We observe that undulations on the sheet grow and de- 
tach to form vortices. At t = 2320, the sheet is ramified 
into a vortex aggregate. Afterwards, sheetlike phase slips 
successively appear forming an expanding self-organized 
superfluid. As well as vortex proliferation, vortex aggre- 
gates sometimes fuse into sheetlike phase slips. These 
temporal changeovers occur on the system scale with 
large-amplitude second sounds. No clear boundary can 
be seen between the self-organized and ordinary super- 
fluid regions, but high-density vortices are undergoing 
large-scale temporal fluctuations in the transition region. 

In Fig. 3, the inhomogeneity of A{r,t) is shown in the 
x-y plane at 2; = 20 taken a,t t = 940. It originates from 
the relation (16) for the entropy and the temperature. 
That is, if ps increases in some regions, the temperature 
there tends to increase if the entropy change is slower. 

In Fig. 4, time evolution of A is shown at three heights 
with X — y — m a longer time range lO'^ < t < 3.4 x lO'^ 
in the same run producing Figs.l and 2. For t < 2000 
the three points are in an ordinary superfluid and are 
nearly uniformly heated. Subsequent zigzag tempera- 
ture changes are due to the passage of platelike defect 
aggregates. A similar plot was given by Moeur et al^, 
where the temperatures at three thermometers attached 
to the side wall became flat during the passage of a self- 
organized region. In our simulation, a train of defects 
eventually arrived at the top and a normal fluid with a 
steep temperature gradient emerged from the top. 
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FIG. 3. (color online). Scaled temperature deviation A{r,t) 
in the x-y plane at z = 20 taken at t = 940 in the same run as 
in Fig.l and Video 1, corresponding to the top plate in Video 
1. A phase slip is being created with considerable temperature 
inhomogeneity, where (A) = -1.22 and {{A - {A)f) = 0.156^ 
with (■ ■ ■) denoting the areal average. 
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FIG. 4. (color online). Scaled temperature deviation A{r,t) 
at z = 100, 140, and 180 with x = y = after application of 
heat flux from the top in the time range 1000 < t < 3400 in 
the same run as in Figs.l and 2. Inset: l^*]^ vs z at t = 2600, 
which is zero at defect points. At this time, self-organized 
and ordinary superfluid regions are separated at 2; ~ 130. 



In Fig. 4, the advancing speed of the self-organized re- 
gion is about 0.1 (in units of ^/i) or 3 cm/s, while it was 
2 X 10~'^cm/s in one of the experimental data (see Fig. 2 
of Ref.l9). Note that the speed of this phase change is 
determined by the rate of heat input to the system. To 
check this aspect, we fixed the bottom heat flux at the 
top value 1 and resumed the integration using the data 
att — 2500 in the run producing Figs. 1-4. Then we could 
realize a dynamical steady state with a self-organized su- 
perfiuid in a lower part and an ordinary superfluid in an 
upper part. Hereafter the times are 2500 plus those after 
the restart. 




-2.7 4.0 
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FIG. 5. Surfaces of I^Pp — 1.2 viewed from two angles in a 
dynamical steady state at t — 3840. Their colors represent 
the temperature deviation A with the color map below the 
plates. A sheetlike phase slip can be seen at the middle, while 
the other defects are aggregates of vortices. 



In Fig. 5, we display snapshots of the surfaces of = 
1.2 at t = 3840 from above and side. The defect structure 
here is composed of a single sheetlike phase slip in the 
middle and platelike vortex aggregates. The vortices are 
highly connected one another forming a network. Some 
vortices start from the phase slip sheet or from the bot- 
tom wall. In this configuration, defects are absent in a 
layer region above the phase slip plane, but above this 
layer a vortex aggregate appears to form an interface be- 
tween the self-organized and ordinary superfiuid regions. 

In Fig. 6, we present the scaled temperature deviation 
A{r,t) vs z with x = y = at three consecutive times 
in the dynamical steady state. It shows considerable dis- 
placive motions of the defect structure on a time scale of 
70. The temperature variations are stepwise at a phase 
slip and are more gradual across vortex aggregates. In 
accord with the balance relation (5), the zig-zag curves 
of A are separated from the fine of r = by 2.5, which 
is our numerical result of Tsos/t. In agreement with this 
result, the experimental formula (6) yields 

r,os/r = 11.2°-*^"/2.5 = 2.85. (23) 



6 




FIG. 6. (color online). Scaled temperature deviation A{r,t) 
vs z (a; = y = 0) at t = 3770, 3840, and 3910 in a dynamical 
steady state in the same run as in Fig. 4. Heat fluxes at the top 
and bottom are fixed at a common value. Significant defect 
movements can be seen. Line of r = is also shown. Inset: 
A vs z at t — 3830, 3840, and 3850 in the ordinary superfluid 
region, indicating standing second sound waves. 



Note that phase slips in the one-dimensional simulations* 
already gave Tgos in agreement vi^ith the experimentii^. In 
the inset, we demonstrate the presence of second sounds 
oscillating on a time scale of 10 (smaller than those in 
Fig.l) in the upper ordinary superfluid region. 

In Fig. 7, we show the phase 9 and the gradient ampli- 
tude I V^'p in the x-z plane at ?/ = using the same data 
as in Fig. 5 at t = 3840. These are cross-sectional pro- 
files illustrating the complex defect structure. Around 
the defects, 9 varies steeply in the upper plate, while 
peaks and ridges appear in the lower plate. Here we 
have |V*p = |^'p|V6'p cx psv"^ around the defects. 



IV. SUMMARY AND REMARKS 

In summary, we have investigated the self-organized 
superfluid state on earth under a fixed downward heat 
fiux Q = 11.2 erg/cm^ s from the top in a cubic cell of 3.2 
mm length. The periodic boundary condition has been 
imposed in the horizontal directions. 

First, at a fixed bottom temperature, we have exam- 
ined how vortices and phase slips appear near the bottom 
and how a self-organized region expands upward. Sec- 
ond, by setting the bottom heat flux equal to the top 
heat fiux, we have realized a dynamical steady state in 
which the self-organized and ordinary superfluids coex- 
ist. The space-time fluctuations are highly enhanced, 
where a sheetlike phase slip often changes into a vortex 
aggregate and vice verse. In fact, in Fig. 5, we can see 
only a single sheetlike phase slip, while numerous vor- 
tex lines form platelike networks. Our numerical value of 
T in Eq.(2) (the distance from the lambda line) in self- 
organized superfluid states is in good agreement with the 




FIG. 7. (color online). Phase 6{r,t) of ^{r,t) (upper plate) 
in the range < ^ < 27r according to the color bar and |V$|^ 
(lower plate) in the x-z plane at j/ = taken at t = 3840, 
where use is made of the same data as in Fig. 5. 



experimental Tgos in Eq.(6). 

We estimate the temperature drop at a phase slip as 

(AT)p = Qip/Mr^os), (24) 

where ip is the thickness of a phase slip (^ ^) and A(tsos) 
is the thermal conductivity in Eq.(9) at r = Tsos- If the 
defects consist of phase slips only, their density np{< 
per unit length satisfies 7T.p(AT)p = T\G in self-organized 
superfiuid states. Use of Eqs.(9), (12), and (24) yields 

npip = GTAA(r,os)/Q = (Tso„/r,os)"^ (25) 

This estimation is consistent with the results in the one- 
dimensional simulations^^'^^. Here the required inequal- 
ity Up < holds for Q > 1. In accord with this, self- 
organized superfiuid states were observed only for Q > 1 
in the experiments^. Notice that the vortex aggregates 
in Fig. 5 also give rise to temperature drops comparable 
to that of a phase slip as in Fig. 6. 

A self-organized superfiuid expands upward if the heat 
input from the top is larger than the heat output from the 
bottom. After its arrival at the top, a normal fiuid begins 
to expand downward (not shown in this paper). There- 
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fore, there can also be coexistence of a self-organized su- 
perfluid and a normal fluid. This coexistence was re- 
alized in our previous simulations^ by cooling a normal 
fluid at the bottom below the transition. In such cool- 
ing, a self-organized superfluid region grows upward to 
reach the top and then an ordinary superfluid region ap- 
pears to expand downward. On heating or cooling from a 
boundary wall without gravity or at large Q, a Hel-Hell 
interface appears and its movement leads to the phase 
transformation between normal fluid and superfluid. 

Phase slips are well-known singular entities observed 
in quasi-one-dimensional superconductor^S, and ^He^. 
However, under an external potential, they can emerge 
as sheets in three-dimensional geometries. The role of 
such a potential becomes relevant on approaching the 
criticality. In the case of gravity, this crossover occurs 
extremely close to the lambda line. Self-organized super- 
fluidity may be realized more easily in an inhomogeneous, 
intense electric or magnetic field, as stated in Sec. I, or 
under high-speed rotation. 

In the literatureii^J^, vortex generation after quench- 
ing of ^He was proposed but it has not been detected 
in subsequent experiments, where crossing of the lambda 
line by a mechanical expansion is supposed. A similar 
quenching method (using the piston effect) was used to 



induce phase separation in fluids near the gas-liquid crit- 
ical point'^— . In a decompression experiment near the 
superfluid transition, however, we should examine how 
the phase transition is induced by sounds. As recently 
reported'^^, first sounds emitted from a boundary wall 
propagate throughout the cell with a stepwise pressure 
decrease at the wave front. They are multiply reflected 
at the boundary walls until dissipated. If the sound speed 
is lO^cm/s and the cell length is 1 cm, the acoustic traver- 
sal time is 10~^s and is much longer than the relaxation 
time of the order parameter except extremely close to 
the lambda line. Second sounds are also emitted from 
the boundary walls with appearance of superfiuid. Fur- 
ther investigation of these processes is thus needed. 
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